Approach to steady-state transport in nanoscale conductors 
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We show, using a tight-binding model and time-dependent density-functional theory, that a quasi- 
steady state current can be established dynamically in a finite nanoscale junction without any inelas- 
tic effects. This is simply due to the geometrical constriction experienced by the electron wavepackets 
as they propagate through the junction. We also show that in this closed non-equilibrium system two 
local electron occupation functions can be defined on each side of the nanojunction which approach 
Fermi distributions with increasing number of atoms in the electrodes. The resultant conductance 
and current-voltage characteristics at quasi-steady state are in agreement with those calculated 
within the static scattering approach. 
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The static scattering approach has been extensively 
used to treat steady-state transport in mesoscopic and 
nanoscopic conductors. The approach, as originally in- 
troduced by Landauer, treats the sample as a scatterer 
between two leads, which are connected adiabatically to 
two infinite electron reservoirs at different local electro- 
chemical potentials i£ The reservoirs are just conceptual 
constructs which enable one to map the non-equilibrium 
transport problem onto a static scattering onei^* 4 How- 
ever, the ensuing steady state may not necessarily be the 
"true" steady state that is reached dynamically when a 
battery discharges across the sample. In addition, the 
static picture says nothing about the dynamical onset 
of steady states, their microscopic nature, or their de- 
pendence on initial conditions. These issues are particu- 
larly relevant in nanoscale structures where some of the 
assumptions of the static approach may hide important 
physical properties pertaining to the true charge dynam- 
ics. 

In this paper we employ an alternative picture of trans- 
port in nanoscale systems in which we abandon the infi- 
nite reservoirs invoked by Landauer. Instead, as recently 
suggested by Di Ventra and Todorov£, we consider the 
current that flows during the discharge of two large but 
finite oppositely-charged electrodes connected by a nano- 
junction. Unlike the static, open boundary approach, 
the present approach permits one to describe the current 
within a microcanonical formalism where both energy 
and particle numbers are conserved quantities. In addi- 
tion, due to the finite and isolated nature of the system, it 
can be demonstrated^ that the total current flowing from 
one electrode to the other can be calculated exactly us- 
ing time-dependent density-functional theory (TDDFT) 6 
provided that one knows the exact functional, regardless 
of whether the system reaches a steady state or not. 

We find that a quasi-steady state current, though last- 
ing only for a limited period of time, can be established 
in the neighborhood of the nanojunction without any dis- 
sipation. This is simply due to the change in the spread 
of momentum of wavepackets as they move into a nano- 
junction and adapt to the given junction geometry. This 
effect occurs roughly in a time At ~ H/AE, where AE is 
the typical energy spacing of lateral modes in the junc- 



tion. For a nanojunction of width w, AE ~ ir 2 h 2 /m e w 2 
and At ~ m c w 2 /ir 2 H. If w = 1 nm, At is of the or- 
der of 1 fs, i.e., orders of magnitude smaller than typical 
electron-electron or electron-phonon scattering times. 7 
We indeed focus on the electron dynamics after the quasi- 
steady state has been established and make a connec- 
tion between this dynamical picture and the Landauer's 
static approach. To this end, we consider a finite three- 
dimensional (3D) model gold nanojunction and a finite 
quasi-one-dimensional (ID) gold wire (see schematics in 
Fig. (IJ. These are the simplest structures for which 
the quantized conductance and current-voltage charac- 
teristics have been computed using the static scattering 
approach 8-9 and have been measured experimentally for 
similar gold quantum point contacts^* Recently, Hors- 
field et al. have shown that a steady current is generated 
in a similar finite atomic chainii Nevertheless the ques- 
tion of whether a steady state can be reached without in- 
cluding any electron-ion interactions remains unanswered 
in their work^ In addition to answering this question, 
we show that one can define two local electron occupa- 
tion functions on each side of the nanojunction. These 
are shifted in energy by an amount which can be inter- 
preted, in the limit of large electrodes, as the "bias" of 
the corresponding open system. These functions depart 
from the equilibrium Fermi distributions by an amount 
which decreases with increasing electrode size. This ver- 
ifies Landauer's hypothesis that "geometrical dilution" 
of wavefunctions is the most important aspect of a reser- 
voir^ However, contrary to previous conclusions^^ we 
show that finite but long one-dimensional leads do not 
need to widen to constitute good "reservoirs" , as long 
as one considers the electron dynamics in the junction 
before the electrons reach the edge of the system. 



We now begin our study by using a simple time- 
dependent tight-binding (TB) model for noninteracting 
electrons where Coulomb interactions and correlation ef- 
fects are absent. Later, we treat the problem using a 
fully self-consistent TDDFT approach in the adiabatic 
local density approximation (ALDA)^£ 
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FIG. 1: Current passing through the junction as a function of 
time for a 3D nanojunction (schematic is shown) calculated 
with a non-interacting TB model. Each of the electrodes con- 
sists of 5 x 5 x 30 atoms arranged in a simple cubic geometry. 
The inset shows the corresponding current for a linear chain 
of N=60, 70, 80 and 90 gold atoms. In both cases E B = 0.2 
eV (see text). 



Consider the iV-site single-orbital TB Hamiltonian 

N N 

fl™ = £ ei |r,)(r j |+t£|r j )<r i+ i|+H.c., (1) 

i— 1 i 

where there is one orbital state \ri) per atomic site 
with energy ti and transfer matrix element t connect- 
ing nearest-neighbor sites>i£ We then prepare the system 
such that half of the system has a deficiency of electrons, 
and the other half has a surplus. This can be done by 
increasing of the sites on one side of the system by 
an energy "barrier" Eb- For the ID wire (see inset of 
Fig. ^| , the interface between the two regions separated 
by the barrier defines the nanoscale "junction" . Taking 
this state as the initial state of the system, we then re- 
move Eb, and let the electrons propagate according to 
the time-dependent Schrodinger equation (TDSE) with 
the time- independent Hamiltonian H TB i& 

Due to the closed and finite nature of the system, the 
total current can be calculated by time differentiating the 
charge accumulated on one side of the system, i.e., 

, N/2 N L 
n— 1 i—l 

Here ip n (t) are the occupied single-electron states that 
are solution of the TDSE, and Nl is the number of sites 
on the left of the junction interface. Summation over spin 
degrees of freedom is implied. 

The onset of a quasi-steady state for a 3D nanojunc- 
tion is shown in Fig. ^ where we plot the current Eq. @ 
as a function of time for Eb = 0.2 eV. In the inset, we 
show that a similar quasi-steady state current develops 



in ID wires of different lengths, where the initial time en- 
ergy barrier forces electrons to change the spread of elec- 
tron momentum, and hence plays a role similar to that 
of the geometric constriction in the 3D case. In all cases, 
the current initially rises rapidly, but quickly settles in 
a quasi-constant value I ss ^ In the ID structures, small 
oscillations are observed which decay in time. The quasi- 
steady state lasts for a time td during which the electron 
waves propagate to the ends of the wire and back. The 
time td is a few femtoseconds for the considered cases, 
and can be made longer by increasing the length of the 
wires (see Fig. ^ inset). We have thus demonstrated nu- 
merically our initial conjecture: in a closed and finite 
nanoscale system, a quasi-steady state current with a 
finite lifetime can develop even in the absence of dissi- 
pative effects. The steady state is a direct consequence 
of the geometrical constriction experienced by the elec- 
tron wavepackets as they propagate through the junction. 
This is in contrast with the conclusion of Ref^, where 
the establishment of a steady state is attributed to a "de- 
phasing mechanism" of the electrons spreading in infinite 
electrodes. 

In order to calculate the conductance of this closed sys- 
tem, we need to define a "bias" . In this non-interacting 
electron problem, the energy barrier Eb seems a natu- 
ral choice. However, that is not completely satisfying as 
it relates to the initial conditions and not to the elec- 
tron dynamics. Let us instead define local occupation 
numbers for electrons in the left and right regions of the 
system. This concept is typically introduced as a start- 
ing point in the static approach to transport. Here we 
would like to define it dynamically. We then project the 
occupation for each eigenstate \Ej) of the Hamiltonian 

H TB , i.e, f(Ej,t) = J2 n K- E j# n (*))| 2 onto thc lcft " and 
right-hand side of the system, 

f{E s ,t) = E| E (E 3 \n)(r t \r(t))\ 2 

n i<N/2 

+E| E {E,\r t )(n\r(t))\ 2 

n i>N/2 

+ E 2Rc { E (E 3 \n)(n\r(t)) 

n ^ i<N/2 

J2 (rmnKnlEA. (3) 

i>N/2 ' 

We denote the first, second, and third term /x,(£^,t), 
fn(Ej,t), and fc(Ej,t), respectively. The quantity 
fc (Ej , t) is the sum of cross terms between the energy 
distribution on the left and on the right region. Note also 
that because the set forms a complete orthonor- 

mal basis, J2j h(Ej,t) = n L , £\ f R (E s ,t) = N - n L , 
and J2j fc{Ej,t) = 0, where Ul is the total charge left to 
the junction at a given time. The quantities fi,(Ej,t) and 
fn(Ej, t), normalized to two electrons per state, are plot- 
ted in Fig. Efl>) for two wires of N = 200 and N = 500 
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FIG. 2: Panel (a) illustrates Afi(t) for a linear chain of N = 
200 (dashed line) and N — 500 (solid line) atoms. Panel (b) 
shows their local occupation functions fL(Ej,t) and fn(Ej,t) 
at a small time t after the onset of current. The inset shows 
the absolute difference between Afi(t — 0) and the initial 
energy barrier Eb = 0.4 eV as a function of N. The dashed 
curve is proportional to the function 1 /N . 



atoms immediately after the onset of the current. 

One might naively think of these functions as broad- 
ened Fermi distributions centered at different "chemical 
potentials", separated by an energy Afi(t). A closer ex- 
amination, however, reveals that they cannot be simply 
fitted to Fermi functions with just an effective thermal 
broadening. Instead, the very functional form of these 
non-equilibrium functions is different from a Fermi distri- 
bution. This is not surprising, as in this finite dynamical 
system electrons spread on each side of the junction and 
are not in any sort of local equilibrium in the electrodes. 
In addition, A/z(t) decreases with time because of the 
transport of electrons from one side of the system to the 
other (see Fig. Uta)). However, as long as we evaluate 
/i(£j, t) and fn(Ej,t) at times much less than ta, these 
functions approach two zero-temperature Fermi distribu- 
tion functions centered at two different energies fj,L and 
and fc(Ej , t) tends to zero for every Ej with increas- 
ing number of sites N in the electrodes. These energies 
can be interpreted as two local "chemical potentials" on 
the left and right side of the system, with A/i = fiR — (ml 
approaching the initial energy barrier Eb with increasing 
N. 

This asymptotic behavior as N — > oo is illustrated 
in Fig. 03Tb) where fi,(Ej,t) and fn(Ej,t) at a small t 
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FIG. 3: Current-voltage (I-V) characteristics at quasi-steady 
state of a finite ID gold wire obtained using the TB approach. 
The corresponding conductance is 1.0 Go- The TDDFT cal- 
culation yields a conductance of 0.99 ± 0.03 Go- 



(0 < t <C id) are plotted for different values of N. The 
inset shows the absolute difference between Afi(t = 0) 
and Eb as a function of N. The difference scales with 
TV in the same way as does the separation of eigen- 
states of H TB close to hl and In three dimensions 
N = N x x N y x N z , and it is then easy to prove, in this 
simple TB model, that Afi(t = 0) approaches (albeit 
"nonvariationally" ) Eb as + jj- + jj- with increasing 
number of atoms in the three different directions (see in- 
set of Fig. [2Jb) for the ID case). It is therefore evident 
that with increasing N, local equilibrium distributions 
can be effectively achieved in the two electrodes without 
inelastic effects: the electron waves moving into these re- 
gions are geometrically "diluted" in a practically infinite 
region of space and therefore do not "disturb" the local 
electron occupation. This is the equivalent of Landauer's 
definition of reservoirs^ Our results, however, show that 
this definition can be extended to one-dimensional elec- 
trodes as well. 

All this discussion allows us to define a conductance 
in this closed system in terms of the current at steady 
state I ss and the value of A/j,(t)/e for N — > oo. The for- 
mer converges very fast with increasing number of atoms, 
whereas the latter is the desired "bias" which, in turn, is 
simply the potential "barrier" Es/e at t = 0. The cur- 
rent I ss as a function of Eb /e is plotted in Fig. |3| The 
corresponding differential conductance is about 1.0 Go 
(Go = 2e 2 /h) at all voltages, in good agreement with 
values obtained from the static approach^ and experi- 
mental observations for similar systems^ 

Finally, we study the onset of quasi-steady states in 
the presence of electron interactions that we describe 
at the mean-field level. We illustrate this point using 
TDDFT within the ALDA2i for ID wires. The corre- 
sponding current is plotted in Fig. 0] for different lengths 
of a finite chain of gold atoms kept at a fixed distance 
of 2.8A aparti^ The lifetime of the quasi-steady state is 
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FIG. 4: Current in the middle of the junction as a function of 
time for a linear chain of N= 20, 40, 60 atoms calculated using 
TDDFT-ALDA. The inset shows the electrostatic potential 
drop along the wire at t — for a wire of 40 atoms. 

short due to the limited system size but clearly increases 
with increasing length of the wire. What is more inter- 
esting, however, is the time for the quasi-steady state to 
set in. The initial transient time is found to be less than 
1 fs, consistent with our original estimate. 

The single-particle Kohn-Sham states&S have no ex- 
plicit physical meaning so that the interpretation of the 
corresponding functions in Eq. is less clear. On the 
other hand, the charge density, and thus, the electrostatic 
potential, are well denned quantities. We therefore define 
the conductance in this closed system in terms of the elec- 
trostatic potential drop between two points inside each 
electrode. 5 As in the case of A/z, the potential drop con- 
verges to the t = value (plotted in the inset of Fig. 
with increasing number of atoms^ The corresponding 
differential conductance is about 0.99 ± 0.03 Go, where 
the average value has been determined from the current 
in the wire with N = 60 atoms at t =lfs. 25 It is worth 
pointing out that when the hopping parameter in the 



tight-binding calculation is chosen to match the drop-off 
time td in the TDDFT calculation for the same num- 
ber of atoms, the initial transient time during which the 
quasi-steady state establishes itself in the tight-binding 
calculation is also less than 1 fs. This observation rein- 
forces the notion that the geometric constriction effect is 
present irrespective of the inclusion of mean-field inter- 
actions. 

Finally, one can obtain an order-of-magnitude estimate 
of the electrode size necessary to observe the quasi-steady 
state. The drop-off time td is roughly given by the time 
it takes for an electron to travel at the Fermi velocity 
of the underlying lattice along the length of the elec- 
trode and back, i.e. td ~ L/vj, where L/2 is the linear 
length of one of the electrodes. Therefore, one should 
observe a quasi-steady state if the time necessary for the 
quasi-steady state to be established At is less than td, i.e. 
L > Vfm c w 2 /ir 2 h. The results illustrated in Figs.nand^] 
are consistent with these crude estimates and show that 
a relatively small number of atoms is necessary to repre- 
sent the electrodes, thus making the present approach a 
practical alternative to standard open-boundary calcula- 
tions of transport. 

In conclusion, we have shown numerically that a quasi- 
steady state can be achieved in a nanoscale system with- 
out dissipative effects, simply owing to the geometrical 
constriction experienced by electron wavepackets as they 
approach the nanojunction. We have also provided a 
practical scheme for dynamical conductance calculations 
in finite nanoscale systems that sheds new light on the 
assumptions of the standard static approach to steady- 
state conduction. The approach is also suited to study 
relatively unexplored effects such as transient phenom- 
ena, time-dependent charge disturbances, uniqueness of 
steady states and their dependence on initial conditions. 
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cial support from the Department of Energy (DE-FG02- 
05ER46204). 
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